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p ABSTRACT 

^ ' Attenuation of high-energy gamma rays by pair-production with UV, op- 

> | tical and IR extragalactic background light (EBL) photons provides a link 

. between the history of galaxy formation and high-energy astrophysics. We 

present results from our latest semi-analytic models (SAMs), which employ 
the main ingredients thought to be important to galaxy formation and evo- 
lution, as well as an improved model for reprocessing of starlight by dust 
to mid- and far-IR wavelengths. These SAMs are based upon a ACDM hi- 
erarchical structural formation scenario, and are successful in reproducing a 
■ large variety of observational constraints such as number counts, luminosity 

| and mass functions, and color bimodality. Our fiducial model is based upon 

a WMAP5 cosmology, and treats dust emission using empirical templates. 
This model predicts a background flux considerably lower than optical and 
near -IR measurements that rely on subtraction of zodiacal and galactic fore- 
grounds, and near the lower bounds set by number counts of resolvable sources 
at a large number of wavelengths. We also show the results of varying cosmo- 
logical parameters and dust attenuation model used in our SAM. For each 
EBL prediction, we show how the optical depth due to electron-positron pair- 
^ | production is affected by redshift and gamma-ray energy, and the effect of 

& . gamma-ray absorption on the spectra of a variety of extragalactic sources. We 

conclude with a discussion of the implications of our work, comparisons to 
other models and key measurements of the EBL and a discussion of how the 
burgeoning science of gamma-ray astronomy will continue to help constrain 
cosmology. The low EBL flux predicted by our fiducial model suggests an 
optimistic future for further studies of distant gamma-ray sources. 
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1 INTRODUCTION 

The extragalactic background light (EBL) is the total- 
ity of light emitted by stars and active galactic nuclei 
* E-mail: rgilmore@physics.ucsc.edu (AGN) over the lifetime of the universe. Today, this per- 
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vasive photon background consists of light emitted at all 
epochs, modified by redshifting and dilution due to the 
expansion of the universe. The bulk of the EBL occurs 
at wavelengths from the near-UV to the far-IR. In the 
UV, optical, and near-IR most of the EBL is due to 
direct starlight, as well as a subdo minant contribution 
from active galactic nuclei (AGN) l|Schirber fc Bullockl 
2003). From the mid-IR to submillimeter wavelengths, 
the EBL consists of reemitted light from dust parti- 
cles, including both continuum thermal radiation and 
line emission fro m polycyclic aromatic hydroca rbons 
(PAH) molecules i|Lagache. Puget fc Dole! 120051 ). The 
background at longer wavelengths is dominated by the 
cosmic microwave background, while shortward of the 
Lyman limit the background flux decreases rapidly due 
to attenuation by neutral hydrogen in stellar atmo- 
spheres and the interstellar and intergalactic media. 

Because the production of the EBL is directly 
linked to the star formation history of the universe, 
limits on the EBL can be used to provide constraints 
on the history of galaxy formation and evolution. Ob- 
servations of the extragalactic sky brightness can con- 
strain the local background, however they do not pro- 
vide information about evolution of the background 
with redshift. Direct sky photometry has been at- 
tempted with a number of instruments, most no- 
tably DIRBE and FI R AS, and also grou nd-based op- 
tical telescopes (e.g. (|Mattila et al.ll201.il )). but these 
type of measurements are subject to considerable un- 
certainties due to th e large foreground so urces that 
must be subtracted (jflauser fc Dwekl 1200 lh . Integra- 
tion over discrete sources seen in galaxy su rveys (e.g. 
iMadau fc Pozzetti|[200ol : iKeenan et alJl2010h is another 
way to estimate the EBL, but one that in principle 
can only provide a lower limit due to the possibility 
of unseen sources beyond the magnitude limits of the 
survey instrument or underestimation of true total lu- 
min osity of galaxies due to light in the faint outskirts 
fe.g. lBernsteinl200^ . Observations with highly sensitive 
satellite instruments have provided us with EBL lower 
limits from galaxy number counts across wide wave- 
length ranges. 

High-energy gamma rays can interact with EBL 
photo ns in electron-positron pair-production i nterac- 
tions (|Nikishovl Il962l ; iGould fc Schrederl Il967l ; Ijellevl 
1966). By effectively removing these gamma rays from 
view, this process has the potential to alter the ob- 
served spectra of extragalactic high-energy sources, and 
increasingly occlude those at higher redshifts. The rapid 
development of ground-based gamma-ray astronomy in 
the past 20 years has l ed to a number of attempts - 
e.g.lDwek fc Krennrichl j|2005l);lAharonian et all (|2006h : 
iMazin fc Raud (120071 ); I Albert et all (|2008l ) - to con- 
strain the EBL based on modification to gamma-ray 
spectra, a method that can provide a measurement of 
the EBL that is independent of direct observation. In 
principle, the cosmological history of the EBL could be 
reconstructed by comparing observations of high-energy 



sources at different redshifts to their intrinsic spectra. 
Unfortunately, the emission mechanisms and intrinsic 
spectra of GeV and TeV sources are still poorly under- 
stood. 

Understanding how the EBL is produced and how 
its spectral energy distribution (SED) evolves in red- 
shift requires an understanding of the sources respon- 
sible for its production. This has been attempted by 
different a uthors using a variety of techniques. As enu- 
merated in lDomuiguez et alJ (|201ll ) (Dll), calculations 
of the EBL fall into four general categories: i) forward 
evolution beginning with initial cosmological conditions, 
such as the semi-analytic models used in this work; ii) 
backwards evolution of the well-constrained present-day 
galaxy emissivity according to some prescription; iii) 
evolution of galaxy properties that are inferred over 
some range in wavelength; iv) direct observation of evo- 
lution in galaxy properties over the redshifts providing 
the major contribution to the background light, a cat- 
egory which describes the empirical method developed 
in Dll. 

The last two of these have become much more 
powerful techniques in recent years due to large- 
scale surveys by ground- and space-based instru- 
ments, especially at UV and IR wavelengths, where 
a great deal of progress has taken place in the 
last decade. Some of the first models to account 
for EBL pr oduction by the evolving ga l axy p opu- 
lation were IMadau. Pozzetti fc Dickinson] (| 1998T ) and 
iFranceschinU (1200111. using HST and I SO data, re- 
spectively, and |Pei. Fall fc Hauserl l|l999l ). who looked 
at chemical enrichment data in Lya systems. A two- 
part paper series by T. Kneiske and collaborators 
( Kneiske. Mannheim fc Hartmannl 2002; K neiske et al.1 
12004 ) computed the EBL and subsequently predic- 
tions for attenuation of gamma-ray sources based 
on a parametrization of the SFR density. These 
models separately include the contribution of the 
LIR G/ULIRG popu la tion. An update to this work 
by iKneiske fc Dole] 

attempts to create a 
model using a similar method that produces a 
minimal back ground, and a similar m ethod was 
employed in Razzaque. Dermer fc Finkel (|2009h and 
iFinke. Razzaque fc Dermei ( 20101 ). The aforementioned 
work of Dll used observed evolution up to z — 4 in the 
K-band luminosity function combined with the evolving 
distribution of 25 galaxy SED types from a multiwave- 
length survey of galaxies to estimate the EBL and its 
evolution. 

Other authors have used backward evolution mod- 
els to predict the EBL. These calculations begin with 
the present day galaxy luminosity function and at- 
tempt to trace this function backwards in time by 
assumin g a functional form for the r edshift evolu- 
tion. In iMalkan fc Steckerl (Il998l . l200ll ). IR luminos- 
ity functions from IRAS were extrapolated backwards 
in redshift using power l aw fu nctions. The model of 
IStecker. Malkan fc Scullvl (|2006h updated this work and 
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computed the EBL below the Lyman limit (13.6 eV) 
fo r two different ca s es of stellar evolution. The model 
of I Rowan- Robinson! (|200lh also utilized a 60 /im evolv- 
ing luminosity function, and a four-component spectral 
model for IR and optical emission. One potential prob- 
lem with this method is that it has difficulty account- 
ing for the emissivity contribution of merger-triggered 
starbursts, believed to contribute an increasing fraction 
of the S FR density and IR emiss ivity with increasing 
redshift. iFranceschini et al . (2001) made an attempt to 
account for this starburst phase in a backwards evolu- 
tion model. IFranceschini. Rodighiero fc VaccarH (2008) 
published a sophisticated model using observed lumi- 
nosity functions, and used it to calculate the EBL and 
gamma-ray attenuation. This model uses evolving lu- 
minosity functions in the near-IR up to z = 1.4 for two 
different galaxy populations (spiral and spheroidal) and 
local luminosity functions for the irregular/starbursting 
population, combined with synthetic SEDs to find the 
total emissivity. 

In forward evolution scenarios such as semi-analytic 
models (SAMs), predictions for the evolution of galaxy 
emissivities are made by beginning from the universe 
in its primordial state and simulating the process of 
galaxy formation. This is considerably more involved 
and challenging than the other methods of estimat- 
ing the EBL, but can provide a degree of insight 
into the fundamental astrophysics processes that de- 
termine the emissivity that is lacking in other ap- 
proaches. Semi-analytic models of structure forma- 
tion based on cold-dark matter (CDM) merger trees 
have been used in several pap e rs by our group to 



predict the EBL. iPrimack et al. J 19991) predicted the 
EBL u sing the SAM described in Somerville fc Primackl 
l|l999l) and ISomerville. Primack fc Faberl |200ll ). Later 



work included improved treatment of absorption and 
reemission of s tarlight by dust and updated cos- 
mological data (IPrimack. Bullock fc Somerville! 120051 : 
IPrimack. Gilmore fc Somervillell2008l ). 

In a co m panion paper to this work, 
Somerville et all (|201ll ). which we will hereafter 
refer to as SGPD11, we present a new semi-analytic 
model based on galaxy formation in a WMAP5 cosmol- 
ogy. This model, which will be summarized in the next 
section, incorporates the physical processes thought 
to be most important in determining the evolution of 
these systems. From the luminosity density calculated 
in this model, we have predicted the evolving EBL out 
to high redshift. In this paper, we address the topic of 
gamma-ray attenuation and show how our estimated 
EBL affects high-energy observations of extragalactic 
sources, and discuss the ability of gamma-ray tele- 
scopes to explore the distant universe. We also present 
a comprehensive comparison to the predictions for 
EBL and gamma-ray opacity that have been proposed 
by a number of recent authors. Our current work is 
an update to IPrimack. Gilmore fc Somervillel (2008), 
which presents an earlier stage of our results using a 



'concordance cosmology' (CACDM), with parameters 
largely consistent with WMAP1. This work is also 
closely related and complementary to iGilmore et alj 
(2009), which used the CACDM model as the basis for 
a prediction of the UV background radiation out to 
high redshift, and therefore emphasized the calculation 
of optical depths for gamma-rays below 200 GeV. 
In that paper we included contributions to the UV 
emissivity from quasars, as well as an account of the 
attenuation of ionizing radiation escaping from galaxies 
and processing by neutral hydrogen in the IGM using 
a radiative transfer code. 

In the following Section, we briefly review the 
key elements of the semi-analytic model presented in 
SGPD11. Results are presented in Section [3] begin- 
ning with a review of key results related to the evolv- 
ing background radiation from SGPD11 in 13.11 In Sec- 
tion 13.21 we show how the population of photons in 
our calculated EBL impacts observations of extragalac- 
tic gamma-ray sources through pair-production interac- 
tions. Section [<01 deals with the comparison between our 
predicted gamma-ray optical depths and observations of 
VHE blazars, and constraints on the EBL that other au- 
thors have derived using high-energy data. In Section [4] 
we compare our EBL model with several others that 
have been proposed using a variety of techniques in re- 
cent years. We conclude in Section [5] with a summary 
of results and a discussion of how current and future 
high-energy observations will continue to constrain the 
EBL. 



2 MODEL 

This section summarizes the semi-analytic model 
that is used to predict the EBL in this work. 
This model is based upo n the models that were 



first presented in Somerville fc Primackl (|l999l ) and 
ISomerville. Primack fc Faberl J200ll). with significant 
new ingredients as described in ISomerville et all l|2008f ) 
(S08). Readers should refer to SGPD11, as well as S08, 
for details. 



2.1 Galaxy formation 

We assume a st andard ACDM universe and a Chabrier 
(|Chabrierll2003l ) stellar initial mass function (IMF) that 
does not evolve in redshift. The model presented in this 
work uses cosmological parameters based on WMAP5, 
including a power spectrum normalization of as ~ 0.82, 
a value that is intermediate between the previous find- 
ings of WMAP1 and WMAP3. This value is within 
la of the rec ently-published value from WMAP7 of 
0.809 ± 0.024 jKomatsu et al.ll201ll ). The SAMs used 
here are based on Monte Carlo realizations of dark mat- 
ter halo merger histories calculated usin g the modified 
Press-Schechter (|Sheth fc Tormenlll999l ) and extended 
Press-Schechter methods. 
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Star formation occurs when gas is accreted by the 
galaxy and becomes available after cooling via atomic 
processes. Feedback from supernovae can heat and eject 
the cold gas within the galaxy. This gas will either 
be deposited in the hot reservoir connected with the 
dark matter halo, or returned to the IGM, depending 
on the wind velocity relative to the virial velocity of 
the halo. Star formation in our model occurs in two 
regimes, quiescent star formation in isolated galaxies 
and merger-driven starbursts. The former is treated us- 
ing a rec ipe based on the empirical Schmidt-Ke nnicutt 
relation (|Kennicuttlll989l ; iKennicutt et al.lll998l '). Merg- 
ers drive gas deep into galactic nuclei, fueling black hole 
growth which power AGN-driven winds. Supermassive 
black holes can also produce radio jets that heat the 
hot halo gas and may eventually shut off cooling and 
eventually lead to a cessation of star formation. This 
"quenching" of star formation tends to occur in mas- 
sive galaxies, which are able to build massive BH, and 
which are accreting gas from hot, tenuous halos rather 
than via cold dense filaments. 

The chemical enrichment and star formation his- 
tory of each galaxy are used to pre dict the total emis- 
sion sp ectrum. We have adopted the lBruzual fc Chariot! 
(2003) stellar population models in this work. 



2.2 Dust extinction and reemission 

Light emitted by stars can be absorbed and reemit- 
ted by dust. In the SAM, dust is modeled as a two- 
component dis tribution, using a mod ified version of the 
prescription of I Chariot fc Falll (|2000l ) . which treats sep- 
arately the dense dust in giant molecular clouds that 
contain newborn stars and the much more diffuse cir- 
rus in the interstellar medium (ISM). Extinction from 
the ISM component is proportional to the density and 
metallicity of cold dust, and a slab geometry is as- 
sumed. Stars younger than f0 7 yr are enshrouded in 
a cloud of dust with optical depth tbc,v = Mbc t~v,o, 
where /ibc = 3 and tv,o is the face-on, V-band ex- 
tinction of the ISM component. To calculate extinc- 
tion at other waveleng ths, we have assum ed a star- 
burst attenuation curve (Calz etti et al . 2000) for the dif- 
fuse dust component and a power-law extinction curve 
A x oc (A/5500 A)", w ith n = 0.7, for the birth clouds 
l|Charlot fc Fal]||2000h . 

We consider two possible normalizations for the ex- 
tinction recipes. In our WMAP5+fixed' model, param- 
eters are constant at all redshifts, and are adjusted to 
match observed relations between UV and IR luminosity 
for nearby galaxies. However, as discussed in SGPDI1, 
we have found that this model has difficulty reproduc- 
ing luminosity functions at higher redshift in the UV- 
and optical bands. Motivated by this finding, we have 
created an 'evolving' model with redshift-dependent pa- 
rameters tuned to match the observed UV- and optical 
luminosity functions at all redshifts where they have 
been measured. In this model, total dust extinction is 



scaled by a factor (1 + z)^ 1 at non-zero redshifts, and 
the opacity and lifetime of molecular clouds is scaled 
by a factor z~ l above redshift 1. Because this evolving 
dust model is found to be more successful at matching 
high-redshift data, we favor this model and will refer to 
it as the 'Fiducial' variant in the next section. 

The reemission of IR light by the dust due to ther- 
mal and PAH emission is estimated in our model using 
templates that describe the spectra of galaxies from the 
mid-IR to submillimeter as a function of the total IR 
luminosity, and are based on observations of galaxies in 
the local universe. Energy absorbed by dust from di- 
rect starlight is redistributed in the infrared according 
to a prescribed SED. These templates are embedded in 
our semi-analytic model, and account for emission at 
wavelengths from a few microns to the sub-millimeter, 
including the emission and absorption lines appearing 
in the PAH region. 

The dust emission templates we hav e used 
here are described in iRieke et al.1 <|2009l . R09), 
and are based on observations of II local LIRGS 
and ULIRGS combined with lower luminosity lo- 
cal systems. A comparison of these template s 



with those of Devriendt. Guiderdoni fc Sadatl (|l999h ; 
Devriendt fc Guiderdonil il2000l), which were used in 
Primack et al.l (1200 it). iPrimack. Bullock fc Somerville] 



( 20051 ). and IPrimack. Gilmore fc Somervillel <|200Sl ). is 
available in SGPD11. Being observationally-based, 
these templates suffer to some extent from starlight con- 
tamination at short wavelengths, and there is also a dis- 
continuity in the galactic SEDs where the templates are 
joined with our stellar synthesis models, at about 4 mi- 
crons. In the wavelength range of 2 to 5 microns we have 
attempted to compensate for these problems by fitting 
a power-law extrapolation to the templates. The overall 
effect of this change to the integrated background light 
that we will present in the next section is minimal. 



3 RESULTS 

3.1 Overview of astrophysical results 

In this section, we present results for the evolving back- 
ground light predicted by our semi-analytic models. In 
addition to our fiducial model, which uses a WMAP5 
cosmology combined with the evolving dust model and 
R09 templates described in the last section, we will also 
show results using employing the 'fixed' dust absorption 
model to facilitate comparison with our older work. Re- 
sults from our CACDM model will also be shown when 
relevant. As in the previous section, much of the mate- 
rial here is a review of key results from SGPD11, and 
readers should refer to this work for further detail. We 
will a lso compare with the work of iDommguez et all 
|201ll ) (Dll), a more observationally-driven model that 
provides a useful contrast to our theoretical approach. 

Predictions for the evolving EBL and gamma-ray 
opacity in the fiducial and fixed models are available 
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online in tabulated forrrQ, and can also be found bun- 
dled with the ArXiv source of this paper. 

3.1.1 Star formation and luminosity density history 

The global star formation rate density arising in our 
WMAP5 and CACDM models is shown in Figure Q] 
along with observational estimates obtained using a va- 
riety of tracers. Our model makes predictions that are in 
agreement with the bulk of data for z < 1, and tend to 
be slightly lower than observed at z ~ 2. All measured 
star-formation rates are subject to significant uncertain- 
ties, as seen in the scatter in results for the plotted data. 
Uncertainties in dust extinction impact all results rely- 
ing on UV luminosity. Measurements of Ha and other 
spectral lines must take into account extinction as well 
as metallicity effects. Other authors have attempted to 
measure star-formation rates based on 24 /im and other 
mid-IR observations of warm dust. These results can be 
affected by AGN contamination, as well as PAH features 
that move in and out of the instrument bandpass with 
changing redshift. All of these uncertainties grow with 
increasing redshift, where our knowledge of dust dis- 
tribution and galaxy SEDs becomes less reliable. The 
star-formation rate density inferred from the UV and 
IR luminosity densities of Dll is also shown here. The 
larger value predicted in this work relative to ours is due 
to the considerably higher far-IR emission in Dll; pre- 
dictions in the optical and near-IR are similar. Above 
redshift 1, star- formation rate density predictions from 
Dll are affected by the assumption made about the 
evolution of different galaxy spectral types. 

The luminosity density in our models is predicted 
by summing over the emission from individual galax- 
ies. Results for total galactic emissivity as a function 
of wavelength and redshift are shown in the right-hand 
panel of Figure [2] In the left panel we show predicted 
luminosity density in the local universe, compared to 
constraints at a number of wavelengths. The local lu- 
minosity density has been extremely well-measured in 
the optical and near-IR by large-scale surveys such as 
SDSS and 2MASS, and this provides a strong constraint 
on any model of the galaxy population. At longer wave- 
lengths, we show how our models fit the local data at IR 
wavelengths, including data from IRAS and SCUBA. 

We have also compared our model with observa- 
tional estimates of the evolving luminosity density at 
a number of different redshifts (Figure [3]). The peak 
emissivity redshift in our model changes depending on 
the wavelength considered. At UV bands, the emission 
closely follows the star-formation rate, which peaks at 
z ~ 2.5 in our WMAP5 model and z ~ 3 in the 
CACDM model. Longer wavelengths include significant 
contributions from progressively more evolved stellar 
populations, and therefore peak at later times. Recent 



1 http : //physics .ucsc . edu/~ joel/EBLdata-Gilmore2012/ 



-0.5 




-2.0 



-2.5 



Figure 1. The star-formation rate density in each of the two 
SAMs over cosmic time. The solid and dotted black curves 
are the predictions from our WMAP5 and CACDM mod- 
els, respectively. The dash-dotted red line is the SFR density 
inferred from the model of Dll, which has been converted 
to a Chabrier IMF. Grey data points are from a compila- 
tion presented in S08; and the dashed line is the estimate of 
iHopkins fc Beacoml (2006). 



evolutionary surveys such as DEEP2 and COMBO-17 
allow us to compare the evolution of galaxy emissiv- 
ity against accurate luminosity density data in several 
bands. Emissivity in the UV (1500 and 2800 A) has 
been seen to increase out to nearly z = 2 l)Dahlen et al.l 
2007). Note the large discrepancy between the evolving 
and fixed dust attenuation models at high redshift in the 
B-band and UV; in the evolving model much less atten- 
uation of starlight occurs in early star-forming galaxies. 
At high redshift our fiducial UV predic tions are some- 
what higher than the measurements of iBouwens et all 
(2007): this is largely due to a substantial contribution 
to the total emissivity from faint galaxies which have 
very little extinctio n in our evolv i ng du st model. 

In the B-band. iDahlen et al.l l|2005h find that emis- 
sion increases out to at least z = 1; this paper makes the 
claim that emissivity in the B- and R-bands is consistent 
with being flat in the interval 1 < z < 2. Results at the 
higher redshifts could be sensitive to the faint end slope 
assumed in calculating the luminosity density. In the K- 
b and we match well the local luminosity measurement 
of iKochanek et all ((2001). Available data at higher red- 
shifts seem to suggest a falloff in emissivity beginning at 
about z ^ 1 which we do not find in our models. As dis- 
cussed in lSomerville et all (l201ll ). our model does seem 
to overpredict the K-band luminosity of galaxies at and 
below L*, beginning at redshifts 1 to 2. A correspond- 
ing overproduction of NIR flux is not seen in the local 
luminosity density or K-band counts. 
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Figure 2. Left: The luminosity density of the local universe. The solid black line is the WMAP5 model, and the dotted line is the 
CAC DM model. Data at a num ber of wavelengths is shown from GALEX (blue circles), SP SS ( red stars;lMontero-Dorta fc Pradal 
|2009|) , 6dF (light blue squares; [jones et alJ|2006h. 2MASS (green star s; ICole et al.ll200ll , and lBell et al.ll2003l) . In the mid- and 
far-IR, the orange squares are from IRAS dSoifer fc Neugebauerlll99ll) , while blue stars are from an analysis of local emissivity 
using data from IRAS, ISO, and SCUBA Ta keuchT^t*aul200lli . Right: Three-dimensional representation of the evolution of 
the luminosity density in our WMAP5 model as a function of wavelength and rcdshift. 



3.1.2 Local EBL flux and discrete sources 

As mentioned in the Introduction, measurements of the 
local (z — 0) EBL generally fall into two categories: di- 
rect sky photometry and integrated counts of galaxies. 
Direct measurements provide an absolute measurement 
of the background light without regard to the sources re- 
sponsible, but require subtraction of foreground sources 
present in the Milky Way and our solar system in order 
to isolate the extragalactic signal. Integration of galaxy 
counts (galaxies per unit sky area at a given magnitude) 
is a way to set firm lower limits on the EBL, although 
the degree to which these measurements converge on the 
true value often remains controversial. The flux from 
faint sources will converge mathematically if the slope 
of the counts plotted on a log number vs flux diagram is 
flatter than unity, or in terms of mag nitudes if a < 0.4 , 
for ln(iV) oc am. As expounded by iBernsteinl l|2007l ). 
photometry of faint galaxies is fraught with difficulty 
in untangling the faint galactic fringes from the back- 
ground, and it is possible to miss 50% or more of the 
light associated with extended sources in simple aper- 
ture photometry. 

Large scale surveys such as the Sloan Digital Sky 
Survey (SDSS), the 6-degree Field survey (6dF) and 
the 2-Micron All Sky Survey (2MASS) have provided 
us with an accurate accounting of the galaxies in the 
local universe, and surveys with the HST have comple- 
mented this data with extremely deep counts. Satellite 
instruments such ISOCAM, IRAC, and MIPS provide 



data in the mid- and far-IR. A detailed presentation 
of galaxy number counts in our models compared with 
data can be found in SGPD11. 

Our prediction for the local EBL is gener- 
ally in agreement with lower limits from inte- 
grated number counts. In the UV, limits from 
iGardner. Brown fe Ferguson! (|2000l ). are considerably 
higher than t he the measurement from GALEX 
jXu et al .1120051 ) . This may be explained by the former's 
use of data from the balloon-based FOCA experiment 
to find bright counts, which were in disagreement with 
those from GALEX at sev eral magnitudes. P reliminary 
Herschel counts data from lBerta et all (|2010h set only a 
weak lower limit on the FIR background peak, and the 
author acknowledges that only about half the total IR 
background is likely being resolved. 

Absolute measurements of the EBL require the re- 
moval of foreground sources, including stars, ISM emis- 
sion, and sunlight reflected from dust in the inner so- 
lar system (often called 'zodiacal' light). The most ro- 
bust direct measurements of the IR background to date 
come from the Diffuse Infrared Background Experi- 
ment (DIRBE) and Far-Infrared Absolute Spectropho- 
tometer (FIRAS) instruments on the Cosmic Back- 
ground Explorer (COBE) satellite, though they are still 
fraught with uncertainty in sk y subtraction (see Fig- 
ure 2 in lHauser fc Dweki lj200ll )). The near-IR flux has 
been calculated from DIRBE observat i ons by a va- 
riety of authors l|Wright fc Reesel |2000| ; IWrighj l200ll ; 
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Figure 3. The luminosity density (integrated luminosity function of sources within a given redshift range) in our models 
vs. redshift at 1500 A, 2800 A, and in the B-, z-, and K- bands (approximately 4500 A, 9130 A, and 2.2 (im, respectively). 
The final panel on the lower right shows the total amount of energy that is absorbed and reradiated by dust at infrared 
energies; units on the y-axis for this panel are solar luminosities per Mpc 3 . The solid black line is our WMAP5 prediction with 
evolving dust, and the dotted line are prior results from our C-ACDM model. Dash-dotted violet shows the predictions from our 
WMAP5 model using fixed dust attenuation parameters. The long -short dashed red lin e is the prediction of Dll. Observational 
data shown here are as follows: 1500 A: Blue sc luares are from iDahlen et al ] ll2007l). red stars are from ISchiminovic h et al.l 
1120051 ). green stars are from iBouwens et al. | ||2007|). an d oran ge circles are from iReddv et alj 1I2OO8I I. The solid purple circle 
is a local measurem ent with GALEX by IWvder et al] l|2005h . 2800 A: Blue squares and the purp le circle are again from 
IDahlen et al.1 d2007t) and IWyder et al. I 11200511 respectively. R ed stars are from iGabasch et al] {2006). B-band: Blue squares 
are from IDahlen et al.l (2005). Deep and Combo-17 data from Faber ct al ] l|2007n " are shown as red star s and open r ed sq uares, 
respectively (these are very similar). This Combo-17 data is an update to t h at ori ginally presented in I Wolf et al] l l2003l) . and 
we show the original points as gree n stars. The work of iMarchesini et al] J2007h is shown as op en purple hexes, z-band: 
Loca l measurements are pr ovided bv lMontero-Dorta fc Pradal j200gT ) fre d) and | B lanton et al] ll2003f) (green). Blue squares are 
from IGabasch et alj l|2006h . K-band : The local determin ation is from iKochanek et al .1 (|200lh ~ High redshift data are from 
iBarro et al] l l2009t) (blue squ ares) and | Ar nputs e t al] l l2007h (open red circ les). Total IR Lum inosity: obser vational estimates 
of the IR emissivity are fr omlCa puti et al] ]2007l) (open blue pentagons) , IReddv et al ] l|2008l l ( green circles'). iRodighiero et all 
l|2010h (purple stars), and lLe Floc'h et al] ]2005l ) (red crosses). 
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Goriian. Wright fc Charvl |2000|; ICambresv et all l200ll : 



Levenson. Wright fc Johnsonl l2007h using foreground 
source subtraction techniques and modeling of the zo- 
diacal light, and has generally yielded high estimates 
in this range compared to number counts. Another no- 
table attempt to measure the near-IR background was 
iLevenson fc Wrightl (|2008l ). which used IRAC data to 
calculate the best-fit flux at 3.6 /im using a profile-fit 
to estimate the light from the unobservable faint fringes 
of galaxies. These results were 70 per cent higher tha n 
those of the aperture method of iFazio et alJ l)2004l ). 
highlighting the large uncertainties that galaxy fringe 
issues can bring to EBL measurement. 

The present-day EBL obtained in each of our mod- 
els is shown in Figure [4] We also show results from Df 1 
for comparison. The local EBL is calculated by inte- 
grating over the luminosity density at all wavelengths 
beginning at z = 7.5, and accounting for the redshifting 
and dilution of photons as the universe expands. The 
EBL at a redshift zq and frequency i/p in proper coor- 
dinates can be written as l|Peeblesl(l993l l 



J(vo,zo) = — 



dl (1 + z ) 3 
dz (1 + z) 3 



e(y, z) dz, (f ) 



where e(y, z) is the galaxy emissivity at redshift z and 
frequency v — vo(l + z)/(l + 20), and dl/dz is the cos- 
mological line element, which is 



1 



dz 



(1 + y/Qm(l + z) 3 +Q A 



(2) 



for a flat ACDM universe. We assume here that the 
EBL photons evolve passively after leaving their source 
galaxies and are not affected by any further interactions 
except for cosmological redshifting. This is an accept- 
able approximation for photons at energies below the 
Rydberg energy of 13.61 eV. Photons above this energy 
are strongly attenuated by neutral hydrogen when leav- 
ing their galaxy of origin. At higher energies, photons 
are also capable of interacting with residual neutral hy- 
drogen and, if sufficiently energetic, neutral and singly- 
ionized helium in the intergalactic medium. The effect 
of these processe s on the ionizing EBL is the topic of our 
previous work in lGilmore et al. r|2009l ). The total flux of 
the integrated EBL, as well as the contributions from 
the optical-near-IR and far-IR peaks and the mid-IR 
valley for each model are shown in Table [1] 



3.1.3 Evolution of the background flux 

A correct determination of gamma-ray opacity at dis- 
tances beyond the very nearby universe, z > 0.05, re- 
quires accounting for the redshift-dependent evolution 
of the background at all wavelengths. The sharply in- 
creasing star formation rate density from z — back to 
z ~ 2, combined with the (z + l) 4 evolution of proper 
flux in redshift means that the background was con- 
siderably more powerful in the recent past, a fact that 



can only be neglected in gamma-ray attenuation calcu- 
lations for the closest extragalactic sources. With ob- 
servations of VHE extragalactic sources now stretching 
out to redshifts of over 0.5, it is important in compar- 
ing different realizations of the EBL that we focus not 
only on the flux at z — but at higher redshifts as well, 
where behavior may be quite different in a given model 
depending on how galaxies evolve. We show how the 
background develops in our models in two ways in Fig- 
ure [5] The top panels show the proper EBL SED from 
different redshifts in the rest frame, for each of our mod- 
els. The bottom panels show the comoving EBL at those 
same redshifts; this is the background that would be 
seen today if all galaxy emissivity had been shut off be- 
low the indicated redshift. ft can be seen in the top plots 
that the EBL photon density was considerably higher in 
the past at all wavelengths. The most striking increases 
from present day levels are in the mid- and far-IR, and 
in the UV. 

Complementary to Figure [SJ in Figure [S] we show 
how the photons populating the IR EBL at various 
wavelengths today were produced as a function of red- 
shift. As expected from our knowledge of obscured star- 
bursting galaxies at high redshift, the mid- and far- 
IR parts of the EBL came into existence considerably 
sooner than the photons that are part of the optical- 
near-IR peak today. Our results are in reasonable agree- 
ment with a recen t survey of submillimeter galaxies 
|Devlin et alJ 120091, which has found that half of the 
background radiation at 250 /im is produced at z > 1.2, 
with this fraction increasing at longer wavelengths. The 
results for the WMAP5 and CACDM models are qual- 
itatively similar, however due to earlier star formation 
in the latter, a greater percentage of photons were in 
place at a given redshift for all wavebands relative to 
the WMAP5 model. Because the measurements shown 
here are unavoidably incomplete at high redshift, the 
fact that our models overpredict the fraction of light in 
place at early times is not necessarily in conflict with 
these results. 

The rapid increase in flux at all wavelengths with 
increasing redshift means that the attenuation per unit 
distance increases a corresponding amount. Therefore, 
gamma rays from more distant blazars suffer more at- 
tenuation than might be expected from the local EBL 
flux. In addition, the SED shape of the EBL changes, so 
a simple ^-dependent scaling factor is not sufficient to 
allow accurate predictions of spectral modification for 
the more distant sources. 



3.2 Gamma-ray attenuation 

The process of photon-photon scattering to electron- 
positron pairs is well understood from quantum elec- 
trodynamics. The basic kinematic requirement for this 
process is that there must be sufficient energy in the 
center-of-mass frame of the two-photon system to cre- 
ate the pair. Including the effect of interaction angle as 
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Figure 4. The predicted z = EBL spectrum from our fiducial WMAP5 model (solid black) and WMAP5+fixed (dash- 
dotted violet) dust parameters, and CACDM (dotted black) models, compared with experimental constraints at a number 
of wavelengths. Dll is shown for comparison in dashcd-dotted red with the shaded area indicating the uncertainty region. 
Data: Upward pointing arrows indicate lower bounds from number counts; other symbols are results from direct detection 
experiments. Note that some points ha ve been shifted slightly in waveleng th for clarity. Lower limits: The blue-violet trian- 
les are results from Hubble and STIS llGardner, Brown fc Fergusoifc oO'), while th e purple open triangles a re from GALEX 
.Xu et alJl2005h. The solid gre en and red triangles are from the Hubble Deep Field jMadau fc Pozzettill2000l l and Ultra Deep 
Field iDolch fc Fergusoiill201lll respectivel y, combined with gr ound based-data, and the solid purple triangle is from a mea- 
surement b y the Large Binocula r Camera {Grazian et aI.l [ 20Qcf ), In the near-IR J, H, and K band s, open violet points are the 
limits from lKeenan et alj ll2010h. Open red tria ngles are from IRAC on Spitzer iFazio et al.ll2004h . and the purple triangle at 
15 (im is fr om ISOCAM jHopwood et al]|201fj) on ISO . The lower l imits f rom MIPS at 24, 7 0, an d 160 fim on Spit zer are 
provided bv lBethermin et al.1 feOlOT) (solid blue) and bv lCharv et al.l J2004T) . iFraver et al.1 J200eft. andlDole et alj ||2006! ') (solid 
gold, open gold, and open green, respectively). Lower limits from Herschel number counts <Berta et alj|2010h are shown as 
solid red triangles. In the submillimeter, limits are presented from the BLAST experiment (green points: [Devlin et alj|20091 ). 
Direct Detection: In the optical, orange hexagons are based on data fro m the Pionee r 10/11 Imaging Photopolarimeter 
iMatsuoka et alj|201 if), which are consistent with the older determination of I Toiler! il983T ). The blue star is a determination 
from lMat^Ua*^t al .1 420111 ). and the triangle at 520 nm is a n upper l imit from the same. The points at 1.25, 2 . 2, and 3.5/^m 
are base d upon DIRBE data with foreground su btraction: [ Wr ightl d2 00ll. dark r ed squares), ICambre sy et al l J200ll. orange 
crosses) , iLevenson fc Wrightl ll2008l. red diamond ) iGoriian. Wright fc Charvl l|2000l . purple open hexes) , IWright fc R eese (2000, 
green square") , and iLe vensonT Wright fc Johnsonl ||20 07|. red as terisks). In the far-IR, direct detection measu rements are shown 
from DIRBE l|Wrightll2004l ; ISchlegeI. Finkbeiner fc Davill99Sl . blue stars and solid red circles), and FIR AS jFixsen et alJll998L 
purple bars). Blue- violet open squares are from IR background measurements with the AKARI satellite l lMatsuura et al1 l2010V 
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Wavelength Range 


WMAP5 (Fiducial) 


WMAP5+Fixed 


CACDM 


Dll 


Optical-near-IR peak (0.1 to 8 /im) 


29.01 


24.34 


26.15 


24.47 


Mid-IR (8 to 50 ftm) 


4.89 


5.16 


5.86 


5.24 


Far-IR peak (50 to 500 /an) 


21.01 


22.94 


24.08 


39.48 


Total (0.1 to 500 /mo) 


54.91 


52.44 


56.09 


69.19 



Table 1. The integrated flux of the local EBL in our models (WMAP5 with evolving and fixed dust parameters, and the 
CACDM model) and the model of Dll. Units are nW/m 2 /sr. 




Figure 5. The history of the EBL in each of our models. The top 2 plots show the background flux at past redshifts in the 
WMAP5 fiducial (left) and CACDM (right) models in standard units. Redshifts shown include z = (solid), z = 0.2 (dotted), 
z = 0.6 (short dashed), z = 1 (long dashed), z = 1.5 (dot-short dashed), 2 = 2 (dot-long dashed) and z = 2.5 (long and 
short dashed); also see the key in the upper- left panel. The bottom two plots show the same quantities, but now evolved to 
present-day (co-moving), allowing easy comparison of the EBL in place at a particular time compared to the total in existence 
at z = 0. 
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Figure 6. Buildup of the present-day background photon 
population in 3 different bands (observed- frame) , showing 
the fraction of the present-day co-moving number density in 
place at a given redshift. Solid black is our fiducial WMAP5 
model, dashed is the WMAP5+fixed, and dash-dotted is the 
m odel of Dll. We have c ompared against the measurements 
of iLe Floc'h et al] | |2009|) (blue) at 24/im and Ijauzac et al.l 
j201lh at 70 (green) and 160 /im (red). For the latter two, 
the extent of the shaded region is the la error bound. 



measured in the cosmological frame, this can be written 

2 



^2E 1 E 2 (1 - cos 6») ^ 2m e 



(3) 



where E\ and E2 are the photon energies and 9 is the 
angle of incidence. We are interested here in cases where 
the target background photon has energies from the far- 
IR ( > 1(T 2 eV) to the Lyman limit ( < 13.6 eV). The 
corresponding gamma-ray energies are therefore in the 
TeV or GeV range. We can rewrite Equation[3]to define 
the minimum threshold energy Eth for a background 
photon to interact with a gamma ray of energy E-,, 

2mlc A 



Eth = -=-7- — - TT". (4) 

E 1 {1 — cos6) 

The cross-section for this proc ess is (|Gould fc Schredeil 
ll967tlMadau fc Phinnevll 19961 ) 

3<jt 



a{E u E 2 ,6 
x 2/3(/3 2 



16 

2) + (3 



(1 



/3 4 )ln 



l-P 



where 



2m? c 4 



£i£ 2 (l - cos( 



(5) 



(6) 



and or is the Thompson scattering cross section. 

The cross section is maximized for center-of-mass 
energies of approximately twice the threshold energy 



2m e c 2 , and falls approximately as inverse energy for 
E 2> Eth- H we also account for 9, we find that the like- 
lihood of absorption is maximized for photons at about 
4 times the absolute threshold energy, with one factor 
of 2 from a and another in going from 9 = n ('head-on' 
configuration) to the most probable angle of interaction 
6 ~ 7r/2. If we assume 9 = 7r/2, then we can define the 
characteristic energy or wavelength for the background 
photons that will most strongly affect a gamma ray of 
energy _E 7 as 



Ebg = 



E', 



= 1.044 



or equivalently, 



X bg = 1.188 



(4) 



fim. 



(7) 



(8) 



Gamma rays at a rest-frame energy above 1 TeV are 
most attenuated by the near- and mid-IR range of the 
EBL, while those in the 200 GeV to 1 TeV regime are 
sensitive to light in the near-IR and optical peak in the 
EBL SED. Below 200 GeV it is mainly UV photons that 
have sufficient energy to cause the pair-production in- 
teraction. Below 19 GeV only background photons with 
energies above the Lyman limit of 912 A have sufficient 
energy to interact at any angle in the rest frame, and 
there is little attenuation due to the paucity of such 
photons (|Obll200ll : [Gilmore et al.ll2009h . 

To calculate the optical depth for a gamma ray ob- 
served at energy U 7 , we perform the integral along the 
line of sight to the target at redshift z, 

1 f z ° dl f 1 
t(Ej,z ) = - dz — / du(l-u)x 

2 Jo dz J-i 

x / dE bg n(E bg ,z)a(E^{l + z),E bg ,e). (9) 

J E mirL 

Where we have 

,-. O A 

Emin = Eth (1+0) = 



E~f(l + z)(l- cost 



to account for the redshifting of the gamma-ray en- 
ergy. Here n(E bg , z) is the proper density of target back- 
ground photons as a function of energy E bg and redshift 
z, and it is shorthand for cos 9. dl/dz is the cosmological 
line element, Equation [5] 

For nearby sources, z 0.05, it is sufficient to use 
the local EBL density n(E bg , z = 0). However, as we saw 
in Figure \5\ both the total power and SED of the EBL 
vary strongly with redshift, and in general it is there- 
fore necessary to understand the evolution of the back- 
ground to correctly compute gamma-ray opacities. The 
rapid increase in flux at all wavelengths with increasing 
redshift to z it 2 means that the attenuation per unit 
proper distance increases a corresponding amount. This 
effect means that gammas from more distant blazars 
suffer more attenuation than might be expected from 
the local EBL flux. In addition, the functional form of 
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Figure 7. The attenuation e~ T of gamma-rays vs. gamma- 
ray energy, for sources at Z = 0.03, 0.1, 0.25, and 0.5, and 
1. Results are compared for our fiducial WMAP5 (solid) and 
WMAP5+fixed (dash-dotted violet) models, as well as the 
model of Dll (red dash-dotted). Increasing distance causes 
absorption features to increase in magnitude and appear at 
lower energies. The plateau seen between 1 and 10 TeV at 
low redshift is a product of the mid-IR valley in the EBL 
spectrum. 




Gamma energy (TeV) 



Figure 8. The gamma ray attenuation edges for the 
WMAP5 (solid black) and WMAP5+fixed (dash-dotted vi- 
olet) models and model of Dll (red dash-dotted). The 
curves show the redshift at which the pair-production optical 
depth t reaches the indicated value for a particular observed 
gamma ray energy. The groups of curves from lower left to 
upper right are the contours for r = 1, 3, and 10. We have 
included thin lines to guide the eye at 50 and 100 GeV. 



the EBL changes, so a simple z-dependent scaling fac- 
tor is not sufficient to allow accurate predictions of spec- 
tral modification for the more distant sources. Using the 
line-of-sight integral, Equation [5] we show in Figure [7| 
the optical depth vs. gamma-ray energy for a variety 
of redshifts. A more general way to show EBL atten- 
uation is to plot the 'attenuation edge' redshift where 
the optical depth reaches a certain value as a function 
of gamma-ray energy, and this is presented out to high 
redshift for 3 values of r in Figure [8] This shows how 
telescopes with lower energy thre sholds will allow us to 
peer deeper into the universe. See lGilmore et all ([2009) 
for a similar plot extending to lower gamma-ray energies 
and higher redshift. 

3.3 Results for TeV blazars 

Today, exploration in the VHE (30 GeV to 30 TeV) 
regime is led by >10m-class imaging atmospheric 
Chere nkov telescopes ( IACTs) in cluding theVERI- 
TAS dMaier et all (20081 ), H.E.S.S. |Hinton|[20o3 ), and 
MAGIC (|Cortinal [20051 experiments, and by the LAT 
instrument on the F ermi gamma-ray sp ace telescope 
jAtwood et all 120091 ) and also AGILE (|Tavani et al.l 
120081 ). 

The Fermi LAT spends most of its time in an-all sky 
survey mode, and with its large area of view is therefore 



an ideal instrument for finding high-energy sources. The 
11-month source catalog lists 685 high-energy source s 
associated with blazar candidates ()Abdo et al.ll2010al ). 
While the Fermi LAT has an energy range of 20 MeV to 
~300 GeV, it has a much smaller effective area than the 
current generation of ground-based instruments, and 
data from the instrument is therefore most useful for 
our purposes at energies below the threshold of these 
IACTs, 50-100 GeV. A detailed analysis of the EBL con- 
straints available from all Fermi observations of blazars 
and GRBs to date was t he subject of a recent paper by 
the Fermi collaboration, lAbdo et al.l (|2010bl ). Current 
limits on the EBL available from Fermi o bservations do 
not co nstrain the UV flux predicted in iGilmore et al.l 
(2009) or in the models presented here. 

In this section and the following section, we will fo- 
cus on the effect of the optical-IR EBL on AGN-type 
sources by IACTs at ^ 100 GeV. Ground-based detec- 
tors searching above 100 GeV have identified 37 extra- 
galactic AGN-like sources at the time of this writing, in- 
cluding 32 BL Lac objects, radio galaxies M87 and Cen- 
taurus A, and the flat-spectrum radio quasars 3C279, 
PKS 1510-08, and PKS 1222+21. With the exception 
of the radio galaxies these objects are all blazars, ac- 
creting AGN which generate tightly beamed relativistic 
jets that are oriented at a small angle relative to our 
line of sight. While they account for the large major- 
ity of detected sources above 100 GeV, BL Lac objects 
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are themselves only a small subset (~20%) of all blazar 
sources, the other 80 percent being flat spectrum radio 
quasars like 3C279. 

3.3.1 Constraints from gamma-ray observations 

While uncertainties and likely variation in the intrinsic 
spectrum of blazars make it impossible to directly link 
the observed spectrum to EBL attenuation, it is possi- 
ble to translate limits on the spectra to EBL constraints. 
The standard assumption in placing limits on the EBL 
from individual spectra is that the reconstructed intrin- 
sic spectrum should not have a spectral index harder 
than 1.5; that is, F > 1.5 where dN/dE oc £T r for 
photon count N, or alternatively dF/dE oc .E _ ( r_1 ) f or 
flux F. This figure comes about both on the basis of ex- 
perimental observations (no observed VHE spectrum is 
harder than this value) and theoretical arguments. The 
standard value for a single-zone synchrotron-self Comp- 
ton (SSC) spectrum is V = (a + l)/2; here —a is the 
spectral index of the shock-accelerated electrons, which 
is not harder than 2.0 in most acce leration models with 
radiative cooling (I Aharonianl 1200 1 ) . This can be inval- 
idated by assuming a non-standard spectrum for the 
electrons; a low energy cutoff in the electron energy will 
lead to inverse-Compton accelerated photons w ith an 
index as low as F = 2/3 (Ka tarzvriski et al1l2006l ). 

The most recent limits on the EBL come from ob- 
servations of blazars at more distant redshifts (z > 
0.1) that have been detected by the current genera- 
tion of ground-based atmospheric Cherenkov telescopes 
(ACTs). Observation by H.E.S.S. of two blazars at 
z=0.165 and 0.186 were used to set limi ts on the near-IR 
EBL based on the F ^ 1.5 criterion (jAharonian et all 
20061); in this case th e maximal limit was the model of 
Primack et al.l (120011 ) multiplied by a factor of 0.45. An- 
other paper by the H.E.S.S. group set constraints from 
blazar 1ES 0229+200 at z=0.1396 (|Aharonian et al.l 
l2007bl ). While this blazar is a closer source than the 
two featured in the 2006 publication, the observed spec- 
trum extended above 10 TeV and therefore probed 
the background in the mid-IR (Equation [8]). In this 
regime, the effect of optical depth on spectral modifi- 
cation is minimal due to the approximate A -1 falloff 
in EBL flux, which leads to a constant photon den- 
sity per logarithmic bin, and therefore an approximately 
constant gamma-ray opacity as a function of energy. 
The most distant source observed at very high ener- 
gies at the time of writing is quasar 3C279 at z=0.536, 
observed by th e MAGIC experime nt during a flare in 
February 2006 jTeshima et al.ll2008l ). The spectrum ob- 
served was quite steep, 4.1 ± 0.7 sta t ± 0.2 sys , and ex- 
tended from about 80 to nearly 500 GeV. An anal - 
ysis of the spectral modification l|Albert et all l2008h 
found that there was little room for an EBL flux in 
the optical higher than one consistent with lower limits 
from nu mber counts, approximately equ i valent to the 
model of IPrimack. Bullock fc Somervillel |2005l ). This 



pape r used a modif i ed ve rsion of the 'best fit' model 
from iKneiske et all (|2004h as the upper limit in the 
optical and near-IR from their finding. An alterna- 
tive analysis of t h e spe ctral deconvolution of 3C279 by 
IStecker fc Scullvl (|2009h disp uted this analysis and ar- 
gued t hat th e higher EBL of lStecker. Malkan fc Scullvl 
(20061 . 120071 ) could still lead to a steep best-fit spec- 
trum. However, this higher EBL is inconsistent with 
Fermi high-redshift blazar observations at the 5cr level 
jAbdo et al.ll2010bl ). 

Another approach to the problem is to attempt 
to constrain the EBL by using spectra from se veral 
sources simultaneously. iDwek fc Krennrichl |2005l ) de- 
rived an upper limit at 60/im by declaring invalid those 
realizations leading to unphysical intrinsic blazar spec- 
tra with sharply rising T eV emissions. More re cently, 
this method was used in iMazin fc Rauel (120071 ). who 
applied constraints from all observed TeV blazars to 
a large number of possible EBL functional forms cre- 
ated using a spline interpolation across a grid in flux 
versus wavelength space. The lower bound of the union 
of excluded models formed an envelope representing the 
highest possible background that does not violate any 
constraints. This was done for 'realistic' and 'extreme' 
bounds of V 1.5 and 2/3 respectively, and provided a 
limit on the EBL from the optical to the far-IR. The lat- 
ter case is motivated by the limiting case of a truncation 
at a low energy bound for the relati vistic electrons re- 
spons ible for the IC component, see iKatarzvhski et al.l 
|2006l ). 

In Figure |9l we show recent upper limits from 
gamma-ray observations in relation to the z — 
EBL from our models. All of our models are gener- 
ally in agreement with these bounds across all wave- 
lengths. While o ur fiducial model do es slightly exceed 
the bounds set bv lAlbert et al.l (|2008l ) in the optical, we 
do not find a conflict with the standard F 1.5 bound 
in Table [2] and we are withi n la agreement with the 
harder spectrum observed bv lAleksic et al.l (|2011bl ). 

It is worth pointing out here that, in general, one 
should use caution concerning these constraints. These 
limits on the present-day EBL do not take into ac- 
count the differences in redshift evolution occurring 
in different EBL models, which becomes increasingly 
problematic for more distant sources. Also, a s men- 
tioned above the limits from lAh aroman et all l|2006D 
and lAlbert et al.l |2008T l assume specific forms for the 
optical peak of the background SED with variable nor- 
malization. The exact normalization of the upper bound 
is dependent upon this choice. The method used by 
Mazin & Raue avoids this second issue, but at a cost 
of more conservative limits resulting from considering a 
finite grid in flux-wavelength space. 

The photon density of the EBL increases with 
wavelength at almost all energies relevant to gamma- 
ray attenuation, and therefore the effect on high en- 
ergy spectra is always a spectral softening. However, 
it is possible that local radiation in the vicinity of 
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Figure 9. Present-day flux predicted in our EBL models, 
compared against upper limits from gamma-ray observa- 
tions. Line types are the same as in Figure [J] Upper limits 
are shown as downward-pointing triangles. Red t riangles are 
observations described in iAharonian et ah and green 

are from MAGIC observations of 3C279 < Albert etalll2008h. 
Other limits shown are the realistic limi ts of lM^zrrTfc^aue! 
2007) (purple), and the analysis of iDwek fc Krennrichl 



20051 ) (orange triangle at 60/^m). The reader should con- 
sult the text for more details and caveats in interpret- 
ing this figure. We also show for comparison many of the 
high estimates of the optical and near-IR EBL f rom direct 
photo metry. The open blue diamonds are from iBernsteinl 
1120071 ) . The points at 1.25, 2.2, and 3.5/^m are base d upon 
DIRBE data with foreground subtraction: IWriehtl d200lf) 
(dark red square s'). iGoriian. Wright fc Charvl l|2000l ) (purple 
open hexes), and lWright fc Rees el teOOCft (green square) . The 
small cy an points are from direct photometry with the IRTS 
satellite I IMatsumoto et alj|2005l) . 



a source could have other effects on the spectrum. 
This is particularly true for FSRQ sources such as 
3C279. As a quasar, 3C279 is a much more powerful 
source at optical and UV wavelengths than BL Lac 
objects. It has therefore been suggested that internal 
absorption from the broad-line region of the quasar 
could harden the spectrum by creating an optical depth 
that decreases with energy over the observed interval 
aronian. Khangulvan fc Costamantd l200ct) , due to 
emissio n in a narrow band o f UV w avelengths. An anal- 
ysis by iTavecchio fc Mazinl ((2009) claimed that while 
significant internal absorption was likely, only the more 
extreme models of the broad line region lead to an ac- 
tual hardening of the intrinsic spectrum, and these mod- 
els lead to a large decrease in flux from absorption, 
by a typical factor of > 10 3 . This effect could poten- 
tially harden intrinsic spectra emerging from AGN be- 
yond the bounds discussed above, but only in limited 



extreme cases. More reasonable models of the local ra- 
diation fields with less total absorption were found to 
leave the spectral index softened or unmodified. 



3.3.2 Spectral modification of known TeV blazars 

We have calculated absorption from each of our EBL 
realizations in the observed spectra of known blazars 
that are approximated by power-law functions, and de- 
termined the approximate power law of the de-absorbed 
spectra. The spectra from these objects are not expected 
to be power laws over large energy ranges. The most 
simple theoretical form of the spectra from SSC emis- 
sion is a double-peaked distribution (when plotted as 
vFj), which arises from synchrotron radiation of lower 
energy photons and inverse Compton upscattering of 
those same photons to gamma-rays. In this model, the 
power law measured at VHE scales is an approximation 
to a section of the inverse Compton peak. 

Also, the effect of gamma-ray attenuation through 
pair production does not in general preserve a power-law 
form, as can be seen in the optical depth plot, Figure 
[7] Quantifying attenuation as a simple modification to 
an intrinsic spectral index is an approximation that is 
only valid when considering short intervals in energy 
and fairly low redshifts. The EBL attenuation has also 
been described as an decaying exponential function in 
energy that affects the spectra above some threshold. 
However, this is a misleading functional description of 
the optical depth. The sharp increase in absorption in 
Figure [7] that appears at multi-TeV energies is caused 
by the rapid increase in photon density as one transi- 
tions from the mid-IR minimum in the EBL SED and 
into the redshift-broadened PAH region and far-IR peak 
(note that our SED is plotted in terms of flux density, 
not number density). This part of the EBL is created 
by re-emitted light from cold dust, much of which orig- 
inates in rapidly star-forming galaxies, and there is no 
reason to believe that this absorption feature would be 
related to an exponential form. The power law and the 
exponential cut-off, which are often used to describe 
gamma-ray spectra, are not amenable to describing the 
full non- linear effects of EBL absorption, which is a line- 
of-sight integral over the evolving photon field. Our op- 
tical depths for nearby sources are relatively straight 
from a couple hundred GeV out to this turnover region, 
so we present results for sources with spectra measured 
in this energy range. 

One other note concerns the integration over bins of 
finite width in energy. As attenuation differs across these 
intervals, it changes the weighting of data and therefore 
the mean within the bin. Properly de-absorbing spec- 
tral data points requires incorporating the optical depth 
into the analysis used to produce the points, and not 
just multiplying by e r at the mean of the bin. Correla- 
tions between the data points must also be accounted 
for in effectively measuring error. The effect of scaling 
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the data with a simple multiplication introduces error 
that is likely to grow with redshift. 

Having warned the reader of these caveats, we 
present results for known blazars seen above ~ 100 GeV 
in Table [2] The results in the table are also shown in 
graphical form in Figure [10] This plot shows the amount 
of change in spectral index after EBL deconvolution 
for a number of blazars, as a function of source red- 
shift. The majority of the objects presented here are 
of the high frequency-peaked BL Lac (HBL) type, with 
the exceptions of intermediate-peaked W Comae and 
3C66A, low-peaked BL Lacertae and S5 0716+714, and 
flat spectrum radio quasars 3C279 and PKS 1222+21. 
Values from this table should only be taken as approxi- 
mate, particularly for more distant sources. To avoid the 
hazards of analyzing data bin-by-bin, as mentioned in 
the previous paragraph, we have based our results here 
on the spectral indices and applicable energy ranges as 
derived by the authors referenced, whenever possible. 
In the cases where the spectrum is claimed to continue 
above the turnover in optical depths seen at several TeV, 
the results become strongly dependent upon the high- 
est energy extent of the fit. When an explicit value is 
not mentioned, we use the highest energy point dis- 
played in the spectrum in the reference. It is argued 
in at least a couple of cases (Markarian 421 and 501, 
l|Konopelko et alj2003h ) that the de-absorbed spectrum 
shows the rollover at the top of the IC peak. Recent 
MAGIC observations of Mrk 501 have detected a spec- 
tral peak at energies which vary in correlation with flar- 
ing activity (| Albert et al.ll2007eh . As simple power-law 
functions do not provide a good fit in this case, we have 
omitted Mrk 501 from our analysis. 

Two spectra in our analysis show unusually hard 
reconstructed spectra. In the case of H 1426+428 

aronian et al.|[2002l ). the reported spectrum that we 
have used does not conform well to a power law, a fact 
that the authors attribute to EBL absorption, and to 
be cautious we have used an upper energy of 6 TeV. 
The sp ectral index of 1ES 0229+200 (jAharonian et all 
l2007bh is sensitive to the highest energies used in the 
calculation, and we assume here an upper energy of 7 
TeV, which is close to the second-highest data point 
presented in the reference. In general, the reconstructed 
spectral indices for spectra extending above a few TeV 
are highly sensitive to the highest energy data point 
considered due to the rapid increase in opacity with en- 
ergies at this scale (Figure 0). 

It is also interesting to compare the effect that dust 
modeling has on spectral reconstruction, by comparing 
the results of the fixed and evolving (fiducial) treat- 
ments in the WMAP5 model. The former produces an 
EBL that is slightly more intense in the mid-IR, while at 
shorter wavelengths the evolving dust treatment of the 
fiducial model allows a higher optical and near-IR inten- 
sity. Higher redshift blazars have generally only been 
seen at lower energies, where gamma-ray attenuation 
is produced by the optical-near IR EBL peak created 
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Figure 10. Here we show the results of Table[2]in graphical 
form. The measured spectral index (r; dN/dE <x E~ r ) and 
redshift of each blazar is shown as a black hexagon with 
error bars, with the index corrected via the Fiducial EBL 
shown as a solid orange point, and that corrected by the 
WMAP5+Fixed model as an open blue point. The horizontal 
dotted line shows T = 1.5, which is typically taken as the 
hardest spectrum possible under usual assumptions. Some 
points have been shifted sideways slightly for readability. 



mostly by redshifted direct starlight. As can be seen in 
Table [2] and Figure 1101 the larger optical and near-IR 
flux of the fiducial model produces more spectral change 
than the fixed model in most blazars. For a few blazars 
that have been seen to multi-TeV energies, the fixed 

model produces a har der slope. 

As mentioned in lAcciari et all (|2009bh . determin- 
ing the redshift of blazars can be difficult due to the 
lack of strong line emission. Blazars S5 0716+714, PC 
1553+113, and 3C66A are cases where VHE spectra 
have been published, but the source redshift remains 
uncertain. With a given background model, gamma-ray 
attenuation can be used to place upper redshift limits on 
these sources (e.g.. IPrandini et al.ll20icil ; lYang fc Wand 
l2010h . Applying the standard T ^ 1.5 constraint, we 
have summarized some findings for these redshift con- 
straints in Table [3] Here, we use the upper la bound 
on source spectral index as the basis for the limit. The 
redshift of 3C 66A was initially determined to be 0.444 
from a single emission line, taken to be magnesium-II, 
and corroborated by a weak Lya detection. Assuming 
this is correct, we find for this blazar a reconstructed 
spectrum that is harder than many others on the list, 
but still significantly softer than the standard limit at 
the la uncertainty bound. The H BL PC 1553+113 has 
been detected by both MAGIC dAlbert et all 1200781 ) 
and H.E.S.S. l|Aharonian et al ] l2008al L but the red- 
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Object ID 


Rcdshift 


Reference 


Experiment 


®L7 




r ob3 


^Fiducial 


Tpixod 


Mrk 421 (+) 


0.030 


Konopclko et al. (2008) 


Whipple 


0.2 


8.0 


2.66 


2.46 


2.46 


Mrk 421 


0.030 


Aharonian et al. (1999) 


HEGRA 


0.5 


7.0 


3.09 ± 0.07 


2.88 


2.88 


1ES 2344+514 


0.044 


Albert et al. (2007d) 


MAGIC 


0.14 


5.4 


2.95 ± 0.12 ± 0.20 


2.71 


2.72 


Mrk 180 ( + ) 


0.045 


Albert et al. (2000b) 


MAGIC 


0.2 


1.5 


3.30 ± 0.70 


3.04 


3.06 


1ES 1959+650 


0.047 


Albert et al. (2006c) 


MAGIC 


0.18 


2.0 


2.72 ± 0.14 


2.46 


2.48 


1ES 1959+650 (-) 


0.047 


Tasliafcrri et al. (2008) 


MAGIC 


0.15 


3.0 


2.58 ± 0.18 


2.34 


2.36 


1ES 1959+650 (+) 


0.047 


Daniel et al. (2005) 


Whipple 


0.38 


18.0 


2.78 ± 0.12 ± 0.21 


1.95 


1.91 


BL Laccrtac 


0.069 


Albert et al. (2007b) 


MAGIC 


0.15 


0.9 


3.60 ± 0.50 


3.25 


3.29 


PKS 0548-322 


0.069 


Supcrina et al. (2008) 


H.E.S.S. 


0.2 


3.0 


2.8 ± 0.33 


2.41 


2.45 


PKS 2005-489 


0.071 


Aharonian et al. (2005a) 


H.E.S.S. 


0.2 


2.5 


4.0 ± 0.4 


3.60 


3.63 


RGB J0152+017 


0.080 


Aharonian et al. (2008b) 


H.E.S.S. 


0.24 


3.8 


2.95 ± 0.41 


2.48 


2.51 


W Comae ( + ) 


0.102 


CoRan (2008) 


VERITAS 


0.15 


2.8 


3.81 ± 0.35 ± 0.34 


3.26 


3.32 


PKS 2155-304 


0.116 


Aharonian et al. (2005b) 


H.E.S.S. 


0.16 


0.70 


3.32 ± 0.06 


2.74 


2.81 


H 1426+428 


0.129 


Aharonian ct al. (2002) 


HEGRA 


1.0 


6.0 


2.60 ± 0.60 ± 0.1 


1.74 


1.67 


1ES 0806+524 


0.138 


Acciari et al. (2009a) 


VERITAS 


0.3 


0.7 


3.6 ± 1.0 ± 0.3 


2.68 


2.79 


1ES 0229+200 


0.139 


Aharonian ct al. (2007b) 


H.E.S.S. 


0.5 


7.0 


2.50 ± 0.19 ± 0.10 


1.45 


1.40 


H 2356-309 


0.165 


Aharonian ct al. (2006) 


H.E.S.S. 


0.16 


1.0 


3.06 ± 0.40 


2.11 


2.23 


1ES 1218+304 


0.182 


Albert ct al. (2000a) 


MAGIC 


0.09 


0.63 


3.00 ± 0.4 


2.28 


2.39 


1ES 1218+304 


0.182 


Fortin (2008) 


VERITAS 


0.16 


1.8 


3.08 ± 0.34 ± 0.20 


2.02 


2.14 


1ES 1101-232 


0.186 


Aharonian et al. (_2006) 


H.E.S.S. 


0.16 


3.3 


2.88 ± 0.17 


1.83 


1.91 


1ES 0347-121 


0.188 


Aharonian ct al. (2007a) 


H.E.S.S. 


0.25 


3.0 


3.10 ± 0.23 ± 0.10 


1.94 


2.03 


1ES 1011+496 (+) 


0.212 


Albert ct al. (2007c) 


MAGIC 


0.12 


0.75 


4.00 ± 0.50 


2.95 


3.10 


S5 0716+714 (+) 


0.31f 


Mazin ct al. (2009) 


MAGIC 


0.2 


0.7 


3.45 ± 0.54 


1.47 


1.74 


PKS 1222+21 (+) 


0.432 


Alcksic ct al. (2011a) 


MAGIC 


0.07 


0.4 


3.75 ± 0.27 ± 0.2 


2.32 


2.60 


3C66A (+) 


0.44f 


Acciari ct al. (2009b) 


VERITAS 


0.2 


0.5 


4.1 ± 0.4 ± 0.6 


1.34 


1.75 


3C279 (+) 


0.536 


Albert ct al. (2008) 


MAGIC 


0.09 


0.48 


4.11 ± 0.68 


2.71 


2.14 


3C279 (+) 


0.536 


Alcksic ct al. (2011b) 


MAGIC 


0.15 


0.35 


3.1 ± 1.1 


0.51 


0.97 



Table 2. Reconstruction of the VHE spectral indices of a number of blazars using our two WMAP5 EBL realizations. r o (, s 
is the index reported by the given reference at energies between and J?g 7 , reported in TeV. These are taken from the 

reference if explicitly stated, otherwise the highest and lowest data points presented are used. In some cases the highest energy 
data point presented has large error bars or does not match well with the power-law fit, and we have opted to use the second 
highest point instead. r Fiducial and r Fixcd are the average intrinsic indices after de-absorption by our two EBL models, over the 
range of energies claimed in the detections. Errors on this quantity are the same as in the observed indices, if provided by the 
author. Plus (+) and minus (-) after the source name are used to signify that the detection was claimed in an abnorm ally high 
or low state; readers should consult the references given for further details. Many of these references were taken from IWagnerl 
l|2008h and also the TeV online catalog, | http ://tevcat. uchicago . edu/\ 



] The redshifts of 3C66A and S5 0716+714 quoted here are considered uncertain; see references for details. 



shift remains unknown at this time. Observations with 
the HST have been unable to find a precise distance, 
but suggest a redshift in the ra nge 0.3 < z < 0.4 
(Tr eves. Falomo fc Uslenghil [2007) . potentially making 
this the most distant VHE-detect ed HBL. Studies of the 
intervening Lya forest with COS l|Danforth et al.ll2010t) 
have more recently found a higher redshift bound, z > 
0.395. Our most stringent constraints with the WMAP5 
model put this bl azar at z < 0.40 with the spectral in- 
dex presented by iPrandini et al.l ((2009), or z < 0.48 if 
the in dex is taken at th e uppe r la level. A similar analy- 
sis bv lMazin fc Goebell ( 2007) using MAGIC data and a 
low-level EBL similar t o the level set by gala xy counts 
found z < 0.69, while [Prandini et al.l d2009h p l aced a 
limit z < 0.67 with the model of iKneiske fc Dolel (|201dh . 
This blazar was also considered in Dll, and an upper 
redshift of z ^ 0.85 ± 0.07 was found. Reconstructing 
the intrinsic spectrum at higher redshifts also was found 
to lead to a break in the power-law shape. Demanding 
that such a break be absent leads to a tighter upper 
limit in this reference, z < 0.42. The weight of this ev- 
idence would seem to suggest a redshift near the lower 
limit set by Danforth. S5 0716+714 is an LBL for which 



a spectrum has been recently reported by MAGIC, and 
was previously detected by EGRET and AGILE. Our 
fiducial bound of z < 0.31 (z < 0. 37 at la) is in agree - 
ment with the range reported by iNilsson et all (|2008h . 
z = 0.31 ± 0.08. 



4 COMPARISON WITH OTHER WORK 

In this section we compare the methodology and 
results of our EBL determination with others in 
the recent literature, including the previous predic- 
tions of earlier versions of our semi-analytic model. 
Our prior prediction for th e EBL , presented in 
IPrimack, Bullock fc Somervillel (2005), used a simi- 
lar semi-analytic model of structure formation to 
that which we have used i n this work and in 
IPrimack. Gilmore fc Somervillel <|200Sl ). The WMAP5 
model presented here has similar normalization to the 
2005 model in the optical and near-IR, leading to low 
flux in these wavebands, with only a small amount 
of light unresolved in the deepest number count sur- 
veys. The differences in the spectral shape of the op- 
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Object ID 


Reference 


Experiment 


E Ll 






z 

Fiducial 


z 

Fixed 


S5 0716+714 (+) 


Mazin ct al. (2009) 


MAGIC 


0.2 


0.7 


3.45 ± 0.54 


<0.37 


<0.42 


PG 1553+113 (+) 


Wagner ct al. (2008) 


MAGIC 


0.09 


0.5 


4.2 ± 0.3 


<0.62 


<0.72 


PG 1553+113 


Prandini ct al. (2009) 


MAGIC 


0.15 


0.6 


4.12 ± 0.17 


<0.48 


<0.51 


PG 1553+113 


Aharonian ct al. (2008a) 


H.E.S.S. 


0.23 


1.3 


4.5 ± 0.32 


<0.48 


<0.52 


3C66A (+) 


Acciari ct al. (2009b) 


VERITAS 


0.2 


0.5 


4.1 ± 0.4 ± 0.6 


<0.53 


<0.56 



Table 3. Here we show the upper redshift limits for 3 blazars with uncertain redshift based on the V ^ 1.5 criterion discussed 
in the last section, and using the upper la uncertainty bound on the spectral indices. Some of these blazars were also shown in 
Table [2] assuming there the most likely source redshift. 



tical peaks are due to changes in the application of 
th e dust absorption prescriptions ; in th is work and 
in IPrimack. Gilmore fc Spmervillel J2008I) we use the 
two-component model of lCharlot fc FallTl^OOOl ). which 
leads to more absorption in the UV and emission in 
the mid- and far-IR. A comparison of the different 
emission templates was presented in SGPD11. The 
CACDM model we have presented features a higher 
level of star-formation, particularly at early times, as 
a result of assuming a larger normalization in the ini- 
tial dark matter power spectrum. In the mid- and far- 
IR, all of our new models produce significantly more 
light due to a larger energy budget from absorbed 
starlight. While the prior prediction was not able to 
match the level of light suggested by the number counts 
with ISOCAM at 15^m or the level of far-IR flux in- 
ferred from DIRBE and FIRAS, all of our new mod- 
els are consistent with constraints from number counts 
in the mid and far infrared and with FIRAS, and are 
near the lower determinations of the DIRBE instru- 
m ent at 140 and 240 fj,m . The AKARI measurements 
of iMatsuura et al.l (12010!) and some DIRB E measure- 
ments jSchlegel. Finkbeiner fc Davisl Il998l ) do suggest 
a larger flux in the far-IR peak. 

The EBL model of 

iFranceschini. Rodighiero fc Vaccarll (120081 ) is based on 
luminosity functions from a variety of survey data. 
Recognizing the need for separate treatments of evolu- 
tion in different wavelength regimes, this model treats 
optical and IR components separately, using the recent 
body of data from Spitzer and other experiments. 
Near-IR luminosity function data up to z = 1.4 is used 
for the spiral and spheroidal populations, while only 
local luminosity functions are used for the irregular 
population. Other local IR data are used to constrain 
other regions of the spectrum. The EBL model pre- 
sented by these authors is similar to ours at most 
wavelengths except the far-IR peak. As their model 
has been derived from the same body of cosmological 
data that our own has been compared against, it is not 
surprising to see similar predictions at low redshift. In 
Figure [11] we show the EBL evolution and gamma-ray 
attenuat ion predictions from our models com p ared with 
those of IFranceschini. Rodighiero fc Vaccaril (|2008[). as 
well a s Dll, model 'C in [Finkc, Razzaqu e fc Dermei 
(2010), both predictions of IStecker. Malkan fc Scully 



(2006*) and the best fit model in iKneiske" et all (|2004D . 
which we will discuss in the following paragraphs. Our 
model is seen to evolve similarly to the Franceschini 
model out to redshift 1. 

The work of IKneiske. Mannheim fc Hartmannl 

(|2002l ) calculated the EBL from the UV to far-IR us- 
ing a 'semi-empirical' method based on measured star 
formation rates and spectral synthesis models. Starlight 
is processed by dust, which is modeled as a blackbody 
with three temperature components. Metallicity is as- 
sumed to increase slowly over cosmic time and an av- 
erage global exti nction curve applied to starlight. A 
follow-up paper, iKneiske et al.l l|2004l ). expanded this 
earlier model into 6 realizations, varying in gas temper- 
ature contribution, star formation rate, and UV escape 
fraction. The 'best-fit' EBL in this paper is consider- 
ably higher than any of our models in both the optical 
and far-IR peaks. The discrepancy in EBL normaliza- 
tion between this model and our own likely originates 
in the star formation rate densities assumed, which 
have a much different functional form than our model. 
Knieske's results are based upon a broken power-law for 
for the star-formation history, with a peak at z=1.2. The 
history predicted by our model is considerably lower 
in this epoch, and does not peak until z » 3 for our 
CACDM model, or z x 2.5 for the WMAP5 model (Fig- 
ure [T}. Thus our models have a lower present-day flux, 
but higher flux at the location of Kneiske's peak and 
at higher redshifts. The use of a blackbody spectrum 
to approximate emission in the PAH region also gives 
their EBL SED a somewhat different shape in the mid- 
IR than we find with our templates which include these 
sharp emission f eatur es. The update to this model in 
IKneiske fc Dole! (|201Cll ) produces a background flux that 
is close to the level seen in discrete number counts from 
the optical to mid-IR, and is similar to our Fiducial 
model. 

The recent m odels presented in 

iFinke. Razzaque fc Derm cr (2010) are based on a simi- 
lar technique to Kneiske's work. These models are built 
from the earlier models in iRazzaque. Dermer fc Finkc 
(2009), in which the contribution to the EBL from 
main-sequence stars was found by computing stellar 
emissivities after assuming forms for the global star- 
forma tion history and IMF. lFinke. Razzaque fc Dermerl 
l|2010l ) expanded these models by including dust rc- 
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A(Angstroms) Gamma energy (TeV) 



Figure 11. Left: Our EBL predictions compared with several recent models from the literature. Solid and dotted black lines 
show the proper flux density from our WMAP5 and CACDM models in the local universe and at z = 1 and z = 2. O ther lines are 
from[Franccschim, Rodighi ero fc Vaccaril | |200S| ) (dashed -dotted blue), the best-fit model o f l lKneiske et aljfeOQ^T ) (long-dashed 
green), Dll (long-short dashed red), and model 'O fromlFinke, Razzaque fc D ernier (20icj) (dashed orange). The baseline and 
fast evolution models of [Stcckcr, Mal kan fc ScuIIvl ||2006| ) are the low and high dotted violet points in the z = panel. Right: 
A comparison of the r 77 = 1 attenuation edges for several m odels. Line types are as in the op posite panel. In this instance, the 
lower set of dotted points denotes the fast evolution model of lStecker. Malkan fc Scul ly ( 2006 ) , while the upper are the baseline 
model. 



absorption and emission, as well as post-main sequence 
stars. The authors f avor t heir model ('C') with a 
iBaldrv fc Glazebrookl J2003I) IMF an d sta r-formation 
histor y from Hopkins fc Beacoml (|2006l ) with lCole et al.l 
(2001) parametrization. The Baldry-Glazebrook IMF 
produces slightly more high-mass stars than the 
Chabrier IMF that we have used. This model has 
a slightly higher normalization than ours in the op- 
tical and near-IR. In the mid- and far-IR, we find 
considerably different SED shapes due to the use of 
different techniques in modeling dust emission, as 
was the case with the models of Kneiske et al.; we 
have used templates which are based on the galaxies' 
total IR emission, while these models assume thermal 
blackbody emission at three fixed temperatures. 

T he models of Stecker and c ollaborators , most re- 
cently IStecker. Malkan fc Scu 113 (|2006l . l2007l ). have ex- 
plored the background using backward evolution mod- 
els. This most recent work proposed two SEDs for the 
EBL, using two different assumptions about the pure 
luminosity evolution of the 60 /im luminosity function. 
The SEDs of all galaxies are assumed to be determined 
by this 60 fim emissivity. The 'baseline' model features 
a pure luminosity evolution multiplier of (1 + z) 3 ' 1 out 
to z=1.4, and constant luminosity from there to z=6. 
The 'fast evolution' model evolves even more quickly, 
as ~ (1 + zf to z=0.8 and ~ (1 + zf for 0.8 < z < 1.5. 



Both of these models are considerably higher than ours 
in the optical and near-IR, with the fast evolution model 
about 50% higher in this range than the baseline; the 
discrepancy in the far-IR with our models is smaller. 
It is difficult to compare our model, which deals with 
galaxies in a system of hierarchically merging dark mat- 
ter halos, with this model, in which it is assumed that 
the local galaxy population just grows brighter with in- 
creasing redshift. Our 60 /im luminosity density is not 
found to increase nearly as quickly as assumed in ei- 
ther of these models; we find that both of our mod- 
els can be well-described by a luminosity density mul- 
tiplier of ~ (1 + z) 1 ' 7 out to z w 1.4 at this wave- 
length. As mentioned in the introduction, the high op- 
tical and near-IR flux of the fast-evolution model puts 
it at odds with t he detection of 3C279 by MAGIC 
(|Albert et al.|[2008 | ). which was disputed by Stecker in 
another analysis ( Stecker fc Scullvi 2009). However, the 
large error on the determined de-absorbed spectral in- 
dex (0.5 ± 1.2), and the possibility o f hardening of 
the spectrum by internal absorption dLiu. Bai fc Mai 
l2008l ; fAharonian. Kh angulvan & Costamantell2008l . but 
see also iTavecchio fc Maz in 2009), make it difficult to 
claim this observation as a firm limit on the EBL. Fur- 
ther observations of this and other high-redshift sources 
will likely improve constraints on flux in the optical EBL 
peak. As mentioned above, both models in this work are 
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in serious conflict with the limi ts set by high-reds hift 
blazars observed by Fermi LAT (jAbdo et al.ll2010bl ). 



A comparison of the model of Dll with our fiducial 
model has been presented and discussed in that work. 
Overall, we find good agreement with Dll from optical 
to mid-IR wavelengths from the local universe to z = 1. 
At longer wavelengths, Dll predicts considerably more 
flux from starlight reprocessed by cold dust in galaxies. 
The large discrepancy in results for the FIR peak in the 
EBL highlights the uncertainties involved in modeling 
this region, an ongoing challenge which we address at 
the end of the next section. Discrepancies in this region 
only have an observable effect on the gamma-ray spec- 
tra for the closest blazars. The star-formation rate im- 
plied in Dll (see Figure 12 in that work) is considerably 
higher than our own, and increases more rapidly with 
redshift between z — and z = 1. At UV wavelengths, 
our model produces more light than predicted in Dll, 
and more absorption of gamma-rays below ~ 300 GeV. 
As shown in SGPD11, the UV luminosity in our model 
is near the largest values measured across redshift in in- 
tegrated luminosity functions, and can be considered a 
maximal prediction at these short wavelengths. 

At nearly all wavelengths we have considered, our 
fiducial EBL SED is near the level of flux resolved 
in discrete backgroun d counts, and we fincl agree- 
ment with the claim of iMadau fe Pozzettil (2000) that 
nearly all light in the optical EBL peak is produced 
by discrete sources. Referring back to Figure 0] we 
recognize two places in our calculated EBL SED in 
which there is tension with observations that do not 
rely strongly on foreground estimates, and which may 
signal shortcomings in our spectral modeling. In the 
UV, we find an EBL lower than calculated using a 
combination of HDF and balloon -based FOCA data 
(IGardner. Brown fc Fergusoij|200Ch . The later GALEX 
experiment, while not capable of surveying to the depth 
of Hubble, found a smaller number of bright counts than 
the FOCA data, likely resul ting from differ ences in cal- 
ibration of the instruments (|Xu et al.ll2005l ). It is there- 
fore possible that the higher Gardner points resulted 
from overestimating the bright counts in their determi- 
nation. Our prediction for the local EBL i s above the in- 
tegrat ed number count me asurements of Grazia net al.l 
1(2008(1 in the U-band and iDolch fc Ferguson! l|201ll) in 
the F435W ban d; in the former th e re is significant dis- 
agreement with lMadau fc Pozz ctti (2000j) at faint mag- 
nitudes, but confirmation that the counts are conver- 
gent. In the near-IR, we find good agre e ment with 
published res u lts fr om IMadau fc Pozzettil l|2000l ) and 
iKeenan et al.l ((2010). Our models fall about la below 
the 5.6 /im lower bound from Spitzer. The fact that the 
5.6 /jm limit is higher than that at 4.5 may cast some 
suspicion on this particular measurement, as there is 
no reason to believe such a spectral feature would ex- 
ist. These lower limits are based upon early 'first-look' 
data, and newer results based on a larger set of survey 
data should soon be available (G. Fazio, private commu- 



nication). Additional sensitivity and survey width may 
be achievable in the 3.6 and 4.5 /jm bands by post- 
cryogenic 'Warm Spitzer' surveys; however the longer 
wavelengt h bands will not be attaina ble at elevated tem- 
peratures |van Dokkum et"alll2007l ). 

Our models lie below the level of direct detection of 
the absolute background by calculations using data from 
Hubble WFPC2, DIRBE, and IRTS (see the discussion 
for Figure Q, however we do find con s istenc y with the 
newer optical results o f [Mattila et al.l (1201 lh . and also 
the Pioneer analysis of iMatsuoka et af l (|201ll l. The low 
significance and large error bars on the HST points of 
iBernsteinl i|2007t l mean that these results should not be 
considered inconsistent with an EBL at the level pro- 
vided by resolved sources. All of our models are at 
least lcr below the flux from any of the near-IR di- 
rect detection calculations we have discussed. Limits 
from gamma-ray observations shown in Figure [5] have 
strongly disf avored the highest levels at th i s rang e. As 
discussed in iLevensoru Wright fc Johnson! i|2007l ). the 
present uncertainty in zodiacal light subtraction may be 
intractable without a new mission to directly study this 
foreground, such as the proposed ZE BRA experiment 
(see discussion in IZemcov et al.ll201lD . While our mod- 
els are consistent within la with number count mea- 
surements by MIPS in the mid- and far-IR, they are 
low compared to the DIRBE measurements of the far- 
IR peak. The zodiacal foreground is a sharply decreasing 
function of wavelength in this regime, and the DIRBE 
points are expected to suffer from less systematic error 
here than in the near-IR, especially at 240 fj,m, where 
our models lies beneath the data. However, our IR tem- 
plates include a somewhat simplified treatment of black- 
body emission at these long wavelengths, and we do not 
expect accurate reproduction in this regime. 

iFardal" et al. 1 (120071 ) compared the possible range of 
EBL flux measurements with observations of the fossil 
mass and star formation rate history of the universe. 
As our semi-analytic model reproduces these 3 observ- 
ables, it is worth discussing our work in the context of 
their claim that a top-heavy or 'paunchy' IMF can best 
fit simultaneously these data. This proposal is based on 
the argument that the low levels of estimated stellar 
mass are difficult to reconcile with the present-day EBL 
flux suggested as being intermediate between the direct 
detections and the lower limits from resolved galaxies, 
which is now kno wn to be in c onflic t with blazar lim- 
its (see Figure |9)| . IFardal et ail (120071 ) created 3 models 
of the EBL based on all available observational limits. 
Their minimal model, with total flux of 50 nW/m 2 /sr, 
is set by resolved number counts and is similar to our 
WMAP5 and CACDM models in the optical and near- 
IR out to the K-band. Their best-fit model, based on a 
compromise between number counts compilations and 
the HST and DIRBE direct detection measurements, is 
substantially higher than our CACDM model. The Ri- 
band number counts are well measured by a number of 
surveys (see Figure 14 in SGPD11) which constrain the 
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amount of stellar mass in the nearby universe. Two fac- 
tors alleviate the discrepancy in our model. Our back- 
ground fluxes are near the lowest levels considered in 
iFardal et all (|2007h . with total fluxes of 56.09 and 54.91 
nW/m 2 /sr for our CACDM and WMAP5 models re- 
spectively (Table [1]), and our global Chabrier IMF pro- 
duces more high-mass stars than the diet-Salpeter con- 
sidered as the standard by these authors. For a near- 
IR flux much higher than our models to not overpro- 
duce the K-band counts, this flux would have to arise 
from a high-redshift population of sources unresolvable 
in our current surveys, which extend to > mag 24. As 
mentioned in the Introduction, there are reasons why 
the star-formation rate measures we compare against at 
high-redshift could be biased high. This interpretation 
favors our fiducial EBL model, which has slightly less 
flux than Fardal's lowest model, and is in fairly good 
agreement with integrated star-formation and observed 
K-counts (see SGPD11). 



5 DISCUSSION 

The EBL presents one of the primary barriers to ex- 
tragalactic gamma-ray astronomy with ground-based 
instruments. Our determination of a fairly low extra- 
galactic background across the optical and near to mid- 
IR, support ed by convergence with alternative meth- 
ods such as iFranceschini. Rodighiero fc VaccarH l|2008f ) 
and Dll as well as new direct measurement techniques 
and recent limits from gamma-ray experiments in the 
optical to mid-IR, is an optimistic prediction for the 
future of the field. The weight of this evidence also in- 
creasingly points to an EBL that is well-known over 
this wavelength range, at a level near that of resolved 
number counts, though many questions remain about 
its redshift evolution. The ability of current- and next- 
generation experiments to detect blazars at larger dis- 
tances is a function of several factors: the luminosity 
function and spectral evolution of these objects, the ef- 
fective area (especially at the lowest energies) and duty 
cycles of these instruments, and the details of the in- 
creasingly uncertain non-local EBL at higher redshifts. 
The field of extragalactic VHE astronomy has grown 
considerably in the last five years, and ongoing progress 
on the instrumentation front suggests that many new 
detections may be coming. 

On Figure [H we have placed vertical lines at 50 and 
100 GeV to facilitate comparison with the gamma-ray 
attenuation edge curves. In the fiducial WMAP5 model, 
the universe does not become optically thick (r > 1) to 
100 GeV gamma rays until z ~ 0.9, and z ~ 1.6 for an 
observed energy of 50 GeV. Our Fiducial model does 
predict slightly higher opaci ties at high redshift com- 
pared to the fiducial model of I Gilmore et al.l (120091 ) ; this 
is due to more UV light escaping high-redshift galax- 
ies in our evolving dust model than in the fixed dust 
model used in the CACDM prediction. Nonetheless, the 



EBL does not become a significant barrier to VHE ob- 
servations at these low energies in either model until 
redshifts considerably higher than those for which AGN 
have presently been detected by ground-based instru- 
ments. At the redshift of the current most distant con- 
firmed source, 3C279, we find r = 0.38 for an observed 
energy of 100 GeV, and r = 1 at 175 GeV. 

As mentioned, the approximation of a local EBL 
in optical depth calculations is only valid for nearby 
extragalactic observations. At redshifts above ~0.3, dif- 
ferences in the evolution of star formation and galaxy 
emissivity begin to have a substantial effect on at- 
tenuation; two different EBL models with the same 
present-day normalization could have widely varying 
behavior at these times. For instance, the Kneiske 
models and Stecker's fast evolution models have star- 
formation history peaks at a lower redshift than our 
models predict. In addition to predicting different re- 
sults for the present-day EBL than our model, the 
evolution with redshift is also quite a bit differ- 
ent in these cases. As no direct observations of the 
EBL are possible at nonzero redshift, predicting at- 
tenuation from sources past these distances must be 
made on the basis of models of galaxy evolution, con- 
strained by surveys of luminosity functions at high red- 
shift. Recent surveys of the non-local universe such 
as DEEP 2 and the multi-wavelength GOODS and 
AEGIS surveys have become powerful tools in con- 
straining the EBL at these distances. Observational 
metho d s such as Franceschini. Rodighier o fc Vaccaril 
(2008j), iFinke. Razzaaue fc Dermerl l)2010l l. and Dll, 
which fall under types (ii), (iii), and (iv) respectively 
in the classification scheme discussed in the Introduc- 
tion, are complementary to our semi-analytic approach, 
which is of the first type. Beyond z ~ 1 to 2, uncer- 
tainties in available star-formation and luminosity data 
become large, and theoretical models will continue to be 
required to understand the production and evolution of 
the background. 

The other impact of the shift to higher redshift ob- 
servations by lower energy-threshold instruments such 
as the H.E.S.S. phase-II upgrade and the future CTA 
experiment is the change in the relevant absorbing pho- 
ton population to UV wavelengths. Our models, and 
the others we reference in the previous section, pre- 
dict a rapid falloff in transmission of gamma rays above 
500 GeV for blazars at the redshift of 3C279. Detecting 
emission at or above 1 TeV from sources at this distance 
will require orders-of-magnitude gains in instrument ef- 
fective area, or observations of flare events with similar 
increases in output. The energy range of primary in- 
terest for these types of sources is going to be 50 to 
500 GeV for the next generation of IACTs, plus lower 
energy data from Fermi. Below 200 GeV, it is the UV 
ba ckground that is prim arily responsible for absorption. 
In I Gilmore et al.l (|2009t ) , we addressed the question of 
this background using four different models varying in 
high redshift star-formation and quasar emissivity; the 
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'fiducial' model in that paper that we consider to be 
the most likely was based on the CACDM model in this 
work. The WMAP5 model with evolving dust parame- 
ters presented in this work is generally consistent with 
this model in its predictions for the evolution of UV 
emission, though it does predict a somewhat stronger 
UV background at high redshift. 

One significant weakness of our present approach is 
our use of templates to describe re-emission by dust at 
mid- and far-IR wavelengths, which is relevant to the 
interpretation of multi-TeV data from nearby blazars. 
This method makes the assumption that galaxies of a 
given bolometric luminosity emit light with a similar 
spectral distribution. As discussed in SGPD11, future 
progress in understanding these wavelengths will likely 
require moving past this assumption. While there is 
much progress to be made modeling this part of the 
background SED, new models of dust will only have a 
substantial effect on our calculation of gamma-ray opac- 
ities for the nearest VHE sources. 
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